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ABSTRACT 

We investigate the topology of the completed 2dF Galaxy Redshift Survey, drawing two 
flux-limited samples of the local Universe from the 2dFGRS catalogue, which contains 
over 220,000 galaxies at a median redshift of z = 0.11. The samples are cut at z = 0.2 
and corrected for selection effects. We then use the three-dimensional genus statistic to 
probe the connectedness of the galaxy distribution on scales ranging from 8 to 20 Mpc, 
and compare these measurements with the analytical result for a Gaussian random 
density field, a generic prediction of inflationary models. We demonstrate consistency 
with inflation on the range of scales considered. We then introduce a parameterisation 
of the analytic genus curve formula that is sensitive to asymmetries in genus number 
as a function of density and use it to demonstrate that such phenomena are ruled out 
with 95% confidence between 8 and 16 Mpc. 
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1 INTRODUCTION 

Despite the success of the ACDM model in explaining the 
observations of cosmological surveys, the goal of understand- 
ing the formation and evolution of the largest structures in 
the Universe remains. Questions about the initial conditions 
of the Universe, from which structure evolves, can be an- 
swered by probing the largest scales with sufficient accuracy 
and density, and using methods of analysis that are ade- 
quately robust. 

Topology has proved to be a valuable test of the prop- 
erties of large-scale structure since the development of the 
genus statistic in t he influential work of J. Richard Gott and 
his collaborators dGott et all Il986l : lHamilton et all Il98rj ; 
iGott et all Il987l . 1 19891 ) and the more recent generalisa- 
tion of this quantity (e.g. iMecke et al.lll994l ; iKerscher et al.l 
2000). The connectedness of the large-scale distribution 
of galaxies is used to probe the initial density field and 
the evolution of the peaks within the field, first as they 
grow linearly with the expansion of the Universe, and later 
non-linearly under gravitational clustering. The topologi- 
cal properties o f this evolution have been investigated i n 
analytical form (jBardeen et alj 1 19861 : IColes fc Jonesl Il99l[ ) , 
where the non-linearit y is frequently t r eated using a per- 
turba tion analysis (e.g. lMatsubaral 1994: Matsub ara fc Sutol 
Il996h . and tested against the observed galaxy distribution 
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with many of the galaxy redshift surveys car ried out to 
date l|Vogelev et al.lll994l : ICanavezes et al.lll998l ). 

The genus statistic has often been used as a diagnos- 
tic to compare the primordial density field with that of 
Gaus sian random phase conditions, a prediction of infla- 
tion l|Bardeen et aflll983l ). To date, no significant deviation 
from Gaussian initial conditions has been detected, though 
a common qualification is the wish for denser sampling over 
a larger volume. 

The most r ecent galaxy redshift surveys, described 
in IColless et ail 1120011) (th e 2dF Galaxy Redshift Survey) 
and lStoughton et all l|2002T) (the Sloan Digital Sky Survey), 
offer an opportunity to surmount this problem. Not only are 
they an order of magnitude larger in volume and in galaxy 
numbers, but the statistical methods used to fairly extract 
information from the raw data have increased in sophisti- 
cation to the point where effectively all observed galaxies 
can be included without bias. The success of these sur- 
veys in analysing such statistical p roperties of the galaxy 
distribution as the p ower spectrum ( Cole et al.l 12005). void 
probability function l|Croton et al.ll2004bl) and~the two-point 
and higher-order cor relation functions ( Hawkins et al.ll2003l ; 
ICroton etaTll2004al ) gives great incentive for a topological 
analysis that will go significantly beyond previous work. 

In this paper, we use the completed 2dF Galaxy Red- 
shift Survey (2dFGRS) to exam ine the topology of the lo- 
cal universe. iHovle et al.l l|2002h have carried out a topo- 
logical analysis of an early release 2dFGRS data set, us- 
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ing the two-dimensional genus statistic (jMelott et alj|l989h 
on a scale of 10 /i _1 Mpc. We wish to extend this analysis 
to take advantage of the full public data release, using the 
three-dimensional genus statistic on scales between 8 and 
20 ft" 1 Mpc. The choice to study a range of scales is an im- 
portant one as it allows for the detection of a systematic 
variation in topology. We use a cosmology with Q m — 0.27, 
Qa = 0.73 and h = 0.72. Hereafter we omit the presence of 
h^ 1 when quoting distances in Mpc. 

Section [2] describes the data set and the preliminary 
analysis it is subject to in order to ensure statistical validity, 
before Section previews the topological genus and explains 
how it is calculated. Section [4] then presents our results and 
the method of error analysis, discussing them within the 
context of the field. 



2 REDSHIFT SURVEY 
SAMPLES 

The Two-degree Field Galaxy Redshift Sur- 
vey l|Colless et al.l [2001) has measured the redshift of 
over 221 000 galaxies extending to a median redshift of 
z = 0.11 and distributed mostly in two wedge-shaped slices. 
These form the independent, contiguous data sets that 
are analysed in this paper. The limiting magnitude of the 
survey in the ftj-band is 19.45 and the survey fields are 
tiled so as to provide a near-uniform sampling rate over the 
relevant regions of the sky. 

The raw data requires k + e correction and the ap- 
plication of masks to account for variations in the lim- 
iting magnitudes and redshift completeness of the survey 
regions. These corrections are performed with the pub- 
lically available software written by Peder Norberg and 
Shaun Cole, based on the ca lculation of the selection func- 
tion in iNorberg et ail |2002l) . Information on the software, 
and the code itself, is provided on the survey website at 
http : //www2 . aao .gov. au/2dFGRS/. 



2.1 Selection function 

As with all redshift surveys, there is a need to ac- 
count for redshift-dependent bias introduced by the fixed 
appa rent magnitude limit . Previous studie s of topol- 
ogy i|Canavezes et all Il998l ; iHovle et al.1 12002| ; iPark et all 
120051 ') have used samples constructed by selecting only 
objects above a fixed absolute magnitude cutoff (i.e. a 
volume-limited sample) , effectively downweighting the over- 
represented population. We have chosen instead to use all 
the galaxies present in the NGP and SGP catalogues, and 
upweight those at high-redshift using the survey selection 
fun ction. 

INorberg et~aH l|2002l ) have calculated the 2dFGRS se- 
lection function based on the fcj-band lu minosity func- 
tion, using the STY (ISandage et al. 19791') and ste pwise 
maximum-likelihood (SWML; lEfstathiou et al.1 Il98cf ) esti- 
mators. The list of galaxy weights generated by the selection 
function is produced by the survey masking software. With 
a weighting value for each galaxy, we constructed purely ap- 
parent magnitude-limited samples from the survey data. 

Figure [1] shows the sample size, sample volume and es- 
timated mean inter-galaxy separation as a function of red- 
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Figure 1. Object counts for both NGP and SGP samples, along 
with the volume of the survey and an estimate of the mean inter- 
galaxy separation. The value of r is calculated incrementally, as 
there will be systematically less points (though they have higher 
weightings) at larger redshift. It is this maximum incremental 
value of the smoothing length that needs to be used in order the 
represent the data fairly. 



shift for both the NGP and SGP slices. This estimate of the 
sampling density is found by computing the inverse number 
density of objects in the survey region — i.e. the volume per 
object — and finding the radius of the sphere about each 
galaxy that would contain this volume. 

1/3 

• > — — = —Kr —■ r = I - — I , (1) 
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N 3 



where N is the total number of objects in the sample, and 
V is the number of survey cells multiplied by the volume 
of each cell (4 3 Mpc). Significantly, there is a break in the 
value of r near z ~ 0.2, corresponding to a drop in the spatial 
completeness of the survey, that represents a natural choice 
of maximum sample size. 



2.2 Smoothing the Samples 

Once the two correctly weighted point distributions have 
been produced, it remains to turn them into a continu- 
ous field amenable to a topological measurement. A natural 
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choice of smoothing kernel is a Gaussian, as it minimises 
the effect of high-frequency 'ringing' that occurs when sharp 
cutoffs are present. Moreover, the form of the genus curve 
of a field sm oothed through a Gau ssian window has been 
calc ulated b y Hamil ton et alj (|l986T ). 

M ar tinez et al, I (|2005l ) have argued that Gaussian 
smoothing is restrictive in the sense that there is fixed range 
of smoothing lengths, outside of which properties of the 
smoothing dominate over the signal. They advocate instead 
a multiscale smoothing process, in particular wavelet anal- 
ysis. We suggest that the selection feature of a fixed-scale 
smoothing process is in fact a desirable characteristic, as the 
comparison of topology on different scales provides informa- 
tion that is important for cosmology. 

The transition between the regimes of linear and non- 
linear structure is smooth and very broad in scale. A charac- 
teristic length occurs where the fractional RMS amplitude of 
density perturbations exceeds unity. This is represented by 
the cosmological parameter a n , the fractional RMS ampli- 
tude between spheres of radius n Mpc, with n = 8 frequently 
used as as is not too far from unity. Indeed, the present 
determination of a s from several independent sources (e.g. 
iTvtler et alj 12004 ) has converged near 0.9. We have used 
the value as ~ 0.89, in line with the clustering analysis of 
this data by 2dFGRS team members. This suggests that the 
characteristic smoothing scale is ~ 10 Mpc. 

An important consideration is that the separation of 
points in the field be well below the smoothing length. The 
value of r, defined in Equation [T] for objects at a given red- 
shift, provides an upper bound on this length. These values 
are shown in Figure [2] The break is illustrated by fitting 
two straight lines to different regions of the plot. The value 
of r for a choice of z max is a good guide to the smoothing 
length that will result in a true field, i.e. not a distribution 
that, while continuous, is really still a series of island peaks 
in an otherwise shallow sea. We decided therefore that as 
it was worthwhile exploring structures on a range of scales, 
selecting 8 Mpc as a lower bound for the smoothing length, 
extending up to only 20 Mpc, in steps of 2 Mpc. To examine 
the topology in this range requires that the choice of z„ lax 
be restricted so as to allow smoothing down to ~ 8-10 Mpc, 

i.e. Zmax ~ 0.2. 

The non-periodic geometries of the 2dF NGP and SGP 
regions permit edge effects in the form of a power loss 
near the survey boundaries, resulting in a systematic un- 
derdensities that have a pronounced effect on topologi- 
cal ineasurements ; _ro_undo this, the method described 
bv lMelott fc Dominikll|l993t) was followed, whereby the edge 
effect is emulated in a separate field of constant density, and 
then ratioed out. To do this, a field with the same geom- 
etry as the survey but of constant positive density (set to 
zero elsewhere) is subject to the same smoothing process as 
the data. The smoothed data field is then ratioed, on a cell- 
wise basis, with the smoothed constant field. Conceptually, 
the edge effect has been factored out, at the cost of adding 
some noise to the topological information near the survey 
boundaries. 

The volume affected by this process can be estimated 
by comparing the power inside the survey region before and 
after smoothing. This is a conservative estimate in the sense 
that even cells with small fractional changes in density, well 
away from the edges of the array, will count toward the 
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Figure 2. Mean inter-galaxy separation r as a function of max- 
imum sample redshift. The break near z ~ 0.2 for both slices is 
the prominent feature. As more galaxies were targeted in the SGP 
slice, the values of r are somewhat lower than for the NGP. 



power transmitted outside the survey region - it is akin to 
measuring the 'equivalent volume' of the smoothing process. 
This transmission increases linearly with smoothing scale; 
for the range of scales we consider, the volume affected for 
the NGP region is 12.7% (smoothing over 8 Mpc) and 26.8% 
(over 20 Mpc), and for the larger SGP slice is 10.0% (8 Mpc) 
and 20.7% (20 Mpc). 



3 TOPOLOGICAL ANALYSIS 

At the heart of the topological analysis is the computa- 
tion of the genu s-threshold density relation, described in 
iGott et"al] (|l989l ). The method calculates the genus number 
of a sequence of two-dimensional surfaces embedded in the 
three-dimensional sample volume. The surfaces are defined 
as the boundaries between the regions above and below a 
fixed critical density, parameterised by the quantity v. Each 
value of v excises a monotonically increasing volume of the 
sample. The fractional volume enclosed by the surface as a 
function of v is 
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dt = -Erf c 



(2) 
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where Erf c is the conjugate error function. This number is 
used to find the critica l density of the d istribution, p c . In the 
CONTOUR3D code jWeinberd [19881 ), a Newton-Raphson 
iteration process is used to determine this value. We have 
instead chosen to compute the critical density by sorting the 
three-dimensional density distibution into a one-dimensional 
array ordered by increasing density, so that, as each element 
of the array represents a standard unit of volume, moving 
n% of the way along the array will give the critical density 
corresponding to enclosing (100 — n)% of the volume. That 
is, p c is the mth element of the sorted arrays, for 



m = (1 - Vf)N, 



(3) 



where N is the number of elements in the arrays (N being 
the same for both). A similar meth od has been employed 
previously bv lCanavezes et all (| 1998T ) . 

The genus number of each two-dimensional surface is 
then measured using the standard CONTOUR3D algorithm, 
with, in particular, the same angle deficit matrix for each of 
the staggered configurations, though we have written our 
own code in Fortran 90 to calculate this. Those cells above 
the critical density are labelled with the value 1; those below 
are labelled 0; those outside the survey region, as defined by 
the ra nge in right ascenion and declination in lColless et al.l 
i|200ll ) are labelled —256. This produced a polygonal sur- 
face on which the Gaussian curvature is compressed into 
the vertices of the polygons. As is conventional, we define g s 
in terms of the integrated Gaussian curvature K, 



g s = g-l = f KdA = — D 



(4) 



for the angle deficit defined by Di = 2ir — £\ Vi, where each 
of th e Vj are the angles a round a vertex. 

IBardeen et alj j 19861 ) have derived the genus number of 
isosurfaces through a Gaussian random field, parameterised 
by the quantity v described in Equation [2] 



g(u;N) = N{l-u 2 ) t 



2 /2 



(5) 



where the normalisation factor N depends on the exact 
choice of the power spectrum of the field. For a featureless 
power law of spectral index n, smo othed by a Gaussian o f 
width A, the normaliation factor is (|Martmez fc Saar]|2002h 



N = 



4tt 2 



3 



3/2 



4tt 2 A 3 



3 + n 



3/2 



(6) 



In practice, the genus curve of the cosmological density 
field has agreed well with this model to first-order. How- 
ever, the curve is sufficiently variable that a pure Gaussian 
random field may be inadequate to give a wholly satisfa- 
tory description. Indeed it is entirely sensible to consider an 
alternative model built upon the foundation of a Gaussian 
random field, but differing in parametric properties that are 
mutually independent and open to physical interpretation. 
Motivated by consideration of the moments of the normal 
distribution, define a 'standardisation' of the genus curve 



».(«/'; P 3 ) = N (1 - u' 2 ) exp (-u' 2 /2) 



where 



(7) 



(8) 



and P3 = (N, fi, a) is a three-parameter vector describing 
the transformation. Modifications from P3 = (N, 0, 1) will 
recentre (/j.) and stretch (a) the genus curve. As a further 
probe of this, we define a fourth parameter a to describe 
the imbalance between the two minima of the genus curve, 
increasing (for positive a) the relative abundance of con- 
nec ted high-density clumps (i.e. a 'meatball' topology, as 
per ISoneira fc Peebles! (|l978f l). 
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Each model is now described by the analogous four- 
parameter vector P4 = (N, p,a,a), and the class of Gaus- 
sian random fields is P4 = (N, 0, 1, 0). 

The determination of the components of P4 proceeds 
by minimising the reduced x 2 function 



1 - 



n — k 



Vi ~ g(vi,N,p,a, a) 



(10) 



where t/; and Oi are the data points and their error bars 
respectively, g(vi; N, p, a, a) is the four-parameter modeQ. 
Here n — k represents the number of degrees-of-freedom (the 
number of independent points less the number of parameters 
in the model). 

For the range of v over which this curve is fit, there 
are 51 data points — as these define surfaces through the 
same density field, they cannot be fully independent, so the 
value of n in the reduced \ 2 ls smaller than the number 
of data points used to fit the model. Moreover, increasing 
this number by measuring over a finer grid of v will not in- 
crease the information contained in the data. It is possible 
to estimate the true number of degrees-of-freedom by mea- 
suring the unreduced x 2 values of the analytic fits (Eq. [5} 
to the genus measurements from realisations of a Gaussian 
random field. The mean value of this sample of the x 2 dis- 
tribution is an unbiased estimator of the number of degrees- 
of-freedom. Using 10 realisations of a 128 3 Gaussian ran- 
dom field, smoothed with a Gaussian kernel of width 8 Mpc, 
and a genus curve with 50 points, the number of degrees-of- 
freedom is 45 ± 7.0 (the la range is generated by bootstrap 
resampling). Each of these fits leaves the amplitude of the 
genus curve as a free parameter — for the four-parameter 
model, we subtract one for each of the remaining parame- 
ters. This defines the effective number of degrees-of-freedom 
to be used in calculating the reduced-x 2 - 



4 RESULTS 

Inflation predicts that the universe will be a Gaussian ran- 
dom field on scales unaffected by gravitationally- induced 
non- linear structure formation (IBardeen et al.lll983l ). If the 
universe is of infinite extent and of finite age, it is a matter 
of looking to large enough scales to see this reflected in the 
galaxy distribution. 

The genus surfaces for the SGP and NGP slices, con- 
structed from flux-limited samples, are shown in Figure [3] 
These two surfaces, which represent large (~ 10 7 Mpc 3 ), in- 
dependent regions of the universe, are in good agreement. 



1 The minimisation is per formed using; the Numerical Recipes 
Fortran 77 amoeba routine llPress et al.lll99"3) . 
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Figure 3. Genus surfaces for the two slices in the 2dFGRS. Both exhibit an apparent asymmetry at small values of A; the degree to 
which this is the case at larger scales is unclear as a result of the decline in amplitude with scale length — an effect that is part of the 
smoothing, rather than representative of a change in the power spectrum with scale. 
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Figure 4. Slices through the genus surfaces for the 2dFGRS. Both show the characteristic asymmetry of a meatball topology on small 
scales, but the trend with scale toward or away from a Gaussian random field is not clear by comparing the curves in this relative fashion. 



They appear to agree with Gaussian random field topol- 
ogy, exhibiting the expected symmetry properties and scale- 
dependence of amplitude. A departure from symmetry is 
apparent on small scales, where the amplitude of the genus 
curves is largest. The asymmetry in the regions of negative 
genus ('troughs') is indicative of a greater number of iso- 
lated high-density regions than isolated low-density regions. 
On larger scales, the asymmetry is not apparent, though in 
the surface representation it is unclear whether this is the re- 
sult of the characteristic reduction in amplitude with scale. 
Figure [4] is a plot of slices through this surface at values of 
A between 8 and 20 Mpc, in steps of 2 Mpc. 

To determine the error bars on our sam ple, w e have 
used the bootstrap resampling technique (jEfronl [l982; 



iBarrow et all Il984l ^. in which we create new samples by 
drawing the full number of objects from the original distri- 
bution, with replacement. The resamples are the same size 
as the 2dF data, but some objects are counted twice or more, 
and some do not appear at all (the probability of drawing 
exactly the data is vanishingly small for a data set of this 
size). The variation between resamples provides a measure 
of the spread in the sampling distribution of the parameter 
we are trying to estimate, i.e. the genus curve of the galaxy 
distribution. 

We generated 40 resamples at smoothing lengths of 8, 
10, 12, 14, 16, 18 and 20 Mpc for both the SGP and NGP 
slices. The standard deviation of the resamples provide error 
bars for each of the curves that were overlaid in Figure [4] 
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Figure 5. Genus curves for the SGP (left) and NGP slices at A = 8 (top), 12 and 16 Mpc. The data points are the crosses with la error 
bars calculated using a bootstrap resampling method. The solid line is the best fit model genus curve for a pure Gaussian random field, 
with the amplitude as a free parameter as described in the body text. The fit is a minimum in reduced \ 2 over the region —3.5 v ^ 3.5. 



The genus curves for 8, 12 and 16 Mpc are shown in Figure[5] 
The best fits of these models are shown in Figures [5] and [6] 

4.1 Comparison with volume-limited samples 

Flux-limited samples receive the full benefit of the catalogue 
sampling, but require an accurately determined selection 
function. A retrospective comparison is possible using sam- 



ples limited by volume, which remove the redshift-dependent 
selection of sources by omitting those too faint to be seen 
throughout the sample v olume (a good descript ion of the 
route to follow is given bv lMartfnez fc Saaij[2002T ). 

The competing effects of increasing survey volume and 
decreasing number density create a peak in the sample size 
as a function of the distance limit - one is integrating pro- 
gressively larger regions of space and smaller regions of the 
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Figure 6. As per Figure[5] where now the model is generated using the four-parameter transformation described in the body text. The 
fit is the minimum in reduced x 2 over the region —3.5 ^ u sC 3.5. 



luminosity function. The choice of redshift cut that max- 
imises the number of galaxies in the sample is z ~ 0.16. 
Figure [7] shows the resulting genus curves for the SGP slice, 
using the same smoothing scales as in Figure [S] The agree- 
ment with those curves is good, though the signal from the 
flux-limited sample is slightly better-defined; the ratio of the 
mean uncertainties (nux-limited:volume-limited) is 1:1.03 at 
8 Mpc, 1:1.15 at 12 Mpc and 1:1.56 at 16 Mpc. Studying the 
progression in error bar size with scale for both sets of curves 



leads to the conclusion that the uncertainties are determined 
by (i) the number of independent subvolumes in each sample 
— fixed by the total volume and the smoothing length — 
as well as (ii) the precision with which the density in each 
subvolume is known — determined by the sampling rate — 
with the former being dominant. At 8 Mpc A is about 1.5 
times the mean separation of points in the volume-limited 
sample (cf. 20 times the flux-limited sample separation), so 
we are just beginning to rub up against this limit. There- 
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name 


effect 





N 


amplitude 


Scales all parts of the genus curve 








equally 


1 


/' 


mean 


Determines position of peak with- 








out altering the form by sliding the 








curve along 


2 


17 


variance 


Broadens or narrows the peaks and 








troughs of the genus curve equally 


3 


a 


'skewness' 


Asymmetrises the depths of the 








troughs 



Table 1. Table of moments of the four-parameter model, listing 
their natural analogues with the normal distribution moments 
and the effect of each of the genus curve. 



fore, the benefit of using flux-limited samples is that one can 
work with a larger volume (the choice of z ma x = 0.2, while 
justified by Figure [2] is by no means exclusive) , as well as 
probe smaller scales. 

Sparser sampling restricts the scales that can be probed, 
and the rapidly reducing galaxy count at redshift above z = 
0.16 removes a significant volume from consideration when 
using samples limited by volume. The flux-limited samples 
are therefore extremely successful at providing information 
over a broad range of scales and throughout a large region of 
the Universe, provided that an accurate selection function 
is available. 

4.2 Measurement of asymmetry 

To assess the significance of the asymmetry, we construct a 
composite probability distribution for the parameters of the 
model by calculating x 2 (Pi) over a grid in four-dimensional 
parameter space. This can be used to determine marginal 
distributions for each parameter, by integrating through 
three of the dimensions. However, it is well worth recog- 
nising the influence of each parameter on the genus curve 
and constructing a physical argument to go with the model, 
though it is only the asymmetric parameters that provide 
information about the relative abundance of high- and low- 
density regions. 

The central moments of the normal distribution are de- 
fined as 

/°° 1 2 2 

( !B -p) m -±=e- ( "-' ,) /2ct dx, (11) 
-oo v27T 

and there is a natural interpretation for the four param- 
eters based around this distribution, which is the part of 
the analytic genus curve function that has been modified 
in our model. Table Q] demonstrates the way in which the 
parameters of the model are analogous to the moments of 
the normal distribution. In fact, the skewness parameter a 
is not really a moment in our definition — properly, a cubic 
term of the form a(x — p,) 3 term is required in the exponent 
to satisfy Equation [11] — but the form we have proposed 
represents the simplest modification available to induce an 
asymmetry. 

As was noted above, the parameters that are most 
relevant here are those that invoke an asymmetry in the 
genus curve, i.e. a and p. We have produced joint proba- 
bility distributions for these parameters by assuming that 



for given (a, fx), the observed data represent independent, 
identicall y distributed realisations, so that, invoking Bayes' 
Theorem l|javnesll2003h . 

p(/i, a\D) oc p(ji, a)p(D\(i, a) oc e -x 2 0' a >/ a (12) 

where p(p, a) represents the prior knowledge that we have 
about the distribution of p and a — as we have none, this 
is set to be uniform. In practice, we are creating the sim- 
plest intuitive probability distribution for these parameters 
given our observations, such that larger values of \ 2 rep- 
resent lower probabilities. To determine this distribution, a 
four- dimensional structure containing the \ 2 value for points 
defined by P is created, and the symmetric parameters are 
marginalised out. Figure [5] shows the distribution on sev- 
eral scales for the model presented in Figure [S] The value 
a = represents a symmetric (skew-free) genus curve, and 
when /i = 0we can interpret this value as corresponding to 
a Gaussian random field that has been squeezed or stretched 
in the wings. Deviations from a = would represent skewed 
genus curves, not consistent with inflation, but the results 
here are consistent with a symmetric initial density field. 
The integrated probability outside the range in these figures 
is negligible — from this the level contours indicating no- 
are inferred. 

Figure [9] shows the dependence of the skewness param- 
eter a on scale for both slices, and for a composite average, 
determined by marginalising through fi. The size of the error 
bars suggest consistency with a — on all scales observed. 
As the scale increases, the smoothing process reduces the 
number of independent points in the sample, and the width 
of the joint probability distribution correspondingly broad- 
ens. Consequently, the strongest constraints on the asymme- 
try of the genus curves occur on shorter scales (8-16 Mpc), 
where jaj < 1.0 and \p\ < 0.8 at the 2a level (i.e. with 
95% confidence). This indicates strongly that high- and low- 
density regions are symmetric at these scales. Above ~ 18 
Mpc, the constraints are weakened considerably, but remain 
indicative of consistency with inflation. 



5 CONCLUSIONS 

We have produced statistically complete flux-limited sam- 
ples of the galaxy distribution using the 2dFGRS and its 
associated selection function. We have measured the topol- 
ogy of the smoothed distribution using the genus statistic 
and compared it to a Gaussian random field. These mea- 
surements are the most precise determination of the genus 
curve for the galaxy distribution to date. We have found 
the data to be consistent with Gaussian random field topol- 
ogy on scales of 8 Mpc and above. We interpret this as a 
measurement of the primordial density field, unaffected by 
gravitational clustering, and so find the galaxy distribution 
to be consistent with the prediction of inflation on scales 
above 8 Mpc. Moreover, we reject the presence of asymme- 
tries in the underlying density field at the 2a level on scales 
between 8 and 16 Mpc. On these scales the topology of the 
galaxy distribution cannot be distinguished from that of a 
Gaussian random field using these methods. 

Future work will focus on observations of the effects of 
non-linear gravitational clustering on topology. This should 
be apparent on scales sufficiently small to have been drawn 
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toward a meatball topology in the time since recombina- 
tion. The genus curve of the smoothed galaxy distribution on 
these scales should exhibit a degree of asymmetry (|a| > 1) 
that cannot be reconciled with an underlying Gaussian ran- 
dom field. The intrinsically small dispersion in genus curve 
measurements associated with short smoothing lengths sug- 
gests that such a result would be made with high precision 

indeed. 

iMatsubaral (|l994T ) has derived an expression for the 
genus curve of a Gaussian random field perturbed by non- 
linear gravitational evolution that is valid in the weakly non- 
linear regime and is parametrised by a, the rms of the den- 
sity co ntrast in the cosmological field. However IMatsubaral 
(2003) demonstrates that for a density field with spectral in- 
dex of n = —1. 5, appropriate to the gal axy distribution on 
below 20 Mpc l|Peacock fc Dodds|[l996l ), the effects of the 
non-linearity on the genus curve are slight. At 8 Mpc, well 
into the non-linear regime, no signi ficant perturbation is ob- 
served in our results, or in those of IPark et alj l|2005l ) using 
the Sloan data. Nonetheless, given data of sufficient preci- 
sion and extending to more non-linear scales, such asymme- 
tries may yet be seen. Understanding and quantifying this 
effect will be an important step in describing the behaviour 
of topological measurements in the non-linear regime. 
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Figure 8. Distribution in probability of the skewness and central tendency parameters for the SGP (left) and NGP slices on scales as in 
Figure[5] The curve is generated by comparing the values of x 2 (Mi °0 after marginalising out the symmetric parameters in the model. On 
all scales, these parameters remain close to zero, consistent with the prediction of inflation, with the constraint on this result weakening 
due to the influx of noise on larger scales. The contours are drawn at 1<t,1.5<t and 2cr. 
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Figure 9. The values of the skewness parameter a for the range of scales we have considered in our analysis, with 2<r error bars, 
calculated from the full- width half- maximum of the distributions of a such as those in Figure [8] The 2dF composite is the average of the 
independent NGP and SGP samples. 
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